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Abstract. We have imaged two bona-fide brown dwarfs with TReCS/GEMINI-S and find mid-infrared excess emission that 
can be explained by optically thick dust disk models. In the case of the young (ai2Myr) Cha Hal we measure fluxes at 10.4yum 
and 12.3^m that are fully consistent with a standard flared disk model and prominent silicate emission. For the a lOMyr old 
brown dwarf 2MASS 1207-3932 located in the TW Hydrae association we find excess emission at 8.7yum and 10.4yum with 
respect to the photosphere, and confirm disk accretion as likely cause of its strong activity. Disks around brown dwarfs likely 
last at least as long as their low-mass stellar counterparts in the T-Tauri phase. Grain growth, dust settling, and evolution of the 
geometry of brown dwarfs disks may appear on a timescale of lOMyr and can be witnessed by observations in the mid-infrared. 
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1. Introduction 

There is mounting evidence that dusty disks accompany brown 
dwarfs (BDs) in their early evolutionary stages similar to 
circum-stellar disks around T Tauri stars (TTS). Based on ex- 
cess emission in the near-infrared (JHK), Muench et al. (2001) 
infer that many substellar objects are initially surrounded by 
circumstellar disks. At longer wavelengths the effects of disk 
emission get more pronounced compared to the (cool) brown 
dwarf photosphere. Jayawardhana et al. (2003) and Liu, Najita 
& Tokunaga (2003) find L' (3.8jum) excess emission in the ma- 
jority of substellar objects in young clusters. The frequency 
of young circumstellar disks appears to be similar in the stel- 
lar and in the sub-stellar regime. Physical properties of dust 
disks around BDs can be better determined with mid-infrared 
(MIR) observations. Exploiting the ISOCAM surveys at 6.7 fim 
and 14.3jum in Chameleon (Persi et al. 2000) and p-Ophiuchus 
(Bontemps et al. 2001), Natta & Testi (2001) and Natta et al. 
(2002) tried to explain the spectral energy distribution (SED) 
of the associated young BD population. They show that the 
MIR excess emission of these objects is generally consistent 
with predictions of simple, passive disk models. This suggests 
far-reaching implications like that the general formation mech- 
anism of stars by core collapse and formation of an accretion 
disk extends into the sub-stellar regime. 

Important physical properties of sub-stellar disks remain 
to be understood. For example, the assumption of hydrostatic 
equilibrium in the vertical disk direction - reasonably accurate 
for many T Tauri disks - implies a flared geometry, and the 
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presence of an optically thin surface layer that is heated by di- 
rect illumination from the central source (Chiang & Goldreich 
1997; Chiang et al. 2001; Natta et al. 2001). In these upper 
layers a prominent emission signature of silicates is expected 
around lOyum, as observed for many disks around Herbig Ae/Be 
(Meeus et al. 2001) and TTS (Meeus et al. 2003; Pryzgodda 
et al. 2003). However, the majority of disk models that actu- 
ally fit SEDs of young BDs apparently require a flat geometry 
(Natta et al. 2002). Alternatively, or in addition to geometri- 
cal effects, also grain growth can weaken, or even extinguish 
the lOyum signature (Bouwman et al. 2001). Unfortunately, the 
poor accuracy and wavelength sampling of the ISO photometry 
hamper to draw unambiguous conclusions. 

New, sensitive, mid-infrared instruments at large telescopes 
allow us to probe the lOjum wavelength range of the SEDs 
of young BDs in much greater detail and with higher accu- 
racy. The first ground-based MIR detection of a disk around 
the young BD candidate ChaHa'2 unequivocally demonstrates 
the absence of a silicate emission feature in its disk (Apai et 
al. 2002). Both, flared and flat disk models, are consistent to 
explain the SED of three young brown dwarfs GY5, GYl 1 and 
GY310 (Mohanty et al. 2004). Apai et al. (2004) present first 
evidence for grain growth and dust settling within the disk of 
the young brown dwarf CFHT BD Tau 4, the best characterized 
BD disk up to date (Pascucci et al. 2003; Klein et al. 2003). 

In this paper we report high-quality MIR photometry of 
two bona-fide BDs, namely ChaHal and 2MASS 1207-3932 
(in the following 2M1207), obtained with the Thermal-Region 
Camera Spectrograph T-ReCS mounted at Gemini-South. A 
disk around ChaHal has already been inferred from ISOCAM 
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Fig. 1. TReCS image of ChaHal at IQAfim. Note the smearing 
due to imperfect tracking. 




Fig. 2. TReCS image of 2M1207-3932 at 10 Afim. Chop and 
nod amplitudes were larger then the field-of view («; 29" x 22") 



observations (Natta & Testi 2001). With the help of our new 
photometric values we refine its disk model. For 2M1207 we 
present the first direct evidence for the existence of circum- 
sub-stellar material around an approximately lOMyr old ob- 
ject. 2M1207 has been identified as BD by optical spectroscopy 
(Gizis 2002), and is related to the TW Hydrae association 
(TWA), a loose group of ~ lOMyr old T Tauri stars at an ap- 
proximate distance of 60 pc to the sun (Webb et al. 1999). The 
lack of a measured K-L' excess did not indicate the presence of 
any circum-sub-stellar material (Jayawardhana et al. 2003), but 
strong emission in Ha and Hel suggested indirectly that this 
object may be a (weak) accretor (Mohanty, Jayawardhana & 
Barrado y Navascues 2003; in the following MJBN). Our new 
measurements - to our knowledge among the most sensitive 
ground-based mid-infrared photometry ever reported - demon- 
strate the persistence of dust around older BDs, on timescales 
at least as long as for low-mass stars. 

2. Observations and Data Reduction 

The imaging observations were performed in service mode 
with TReCS (Telesco et al. 1998) mounted at Gemini-South 
on Jan. 2, 2004 (ChaHal) and on Jan. 27, 2004 (2MASS1207- 
3932) under clear weather conditions. Its pixel scale of 0.09" 
results in a field of view of ^ 29"x22". The same throw ampli- 
tudes, but opposite directions, were used for chopping and nod- 
ding. The observations on Jan. 2 were obtained with throw am- 
plitudes of 10 arcsec, and the positive and two negative beams 
are all located in the field of view. The source appears only in 
the final, co-added frame, because slight tracking eiTors of the 
telescope introduced some image smearing of about 2-3 pix- 
els. A more accurate correction of the tracking errors was not 
possible due to the faintness of ChaHal, and a "shift and add" 
co-add of individual chopping pairs was not possible (in con- 
trast to the case of BD-Tau4 discussed in Apai et al. 2004). 
Fig- n shows the final frame with the detection. 2M1207 was 
observed with 15 arcsec throw amplitudes, and only the overlay 
of the two positive beams can be extracted in the final image. 
No tracking errors were evident, and the target appears point- 
like in the final image (Fig. ^. A summary of the observing 
log, including total integration times, is given in Tabled] 



Table 1. Log of the TReCS observations. 



RJD 


Object 


Filter 


Airmass 


Int[s] 


Calibrator 


53006.74 


ChaHal 


Sil-7.7 


1.6 


782 


HD92305 


53006.70 


ChaHal 


Si4-104 


1.6 


912 


HD92305 


53006.79 


ChaHal 


Si6-12.3 


1.5 


1216 


HD92305 


53031.84 


2M1207 


Sil-7.7 


1.0 


304 


HDl 10458 


53031.85 


2M1207 


Si2-8.7 


1.0 


608 


HDl 10458 


53031.81 


2M1207 


Si4-104 


1.0 


782 


HDl 10458 



Final co-added images were obtained from TReCS pipeline 
processed images and from manual merging of individual 
frames, with consistent results. Source count-rates were ex- 
tracted with standard aperture photometry. Curve-of-growth 
methods were applied to find the optimal aperture radius that 
maximized the signal-to-noise ratio of the extracted source. 
That radius was then applied for the calibration source to de- 
termine the count-rate to flux conversion factor. Error estimates 
are based on the formal error of the background noise weighted 
over the source extraction region. When positive and negative 
beams were available, the fluxes were found to be consistent 
with these formal error estimates. The stellar parameters as- 
sumed, and mid-infrared fluxes observed are summarized in 
TableEl 

Given the effective temperatures and bolometric luminosi- 
ties indicated in Table |2 the purely photospheric contribution 
to the 10.4yum flux density is about ImJy for ChaHal and 
1 .5mJy for 2M1207. This estimate is based on BD atmospheric 
models from the Lyon group (e.g. Allard et al., 2000, 2001). For 
the following, we use their published tables that include full 
dust treatment (both condensation and opacities) and use the 
TiO and H2O line lists from Nasa AMES. We refer to a surface 
gravity of log(^) = 3.5 and solar metallicity. The atmospheric 
fluxes in the N-band in these models do not vary by more then 
1 % from a blackbody flux assuming the same Tg/f (in contrast 
to the near-IR, where the deviations can be huge). Hence, both 
our measurements indicate significant excess that cannot come 
from a BD photosphere alone. An excess criterium based only 
on NIR photometry (like E(K - L')o used e.g. in Jayawardhana 
et al. 2003) has very limited value to signal the presence of 
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Table 2. Stellar parameters and mid-infrared fluxes, l.ljjm upper limits, and 8.7/im, 10.4yum, 12.3/im measurements are from 
this work. 



Object 


SpT 


[K] 


Lira; 
[Lo] 


M 

[Mo] 


age 
[Myr] 


[mJy] 


G.lfiin 
[mJy] 


1.1 lim'' 
[mJy] 


S.lfim 
[mJy] 


10.4jum 
[mJy] 


U.Sfim 
[mJy] 


I4.3ijm 
[mJy] 


ChaHal 

2M1207 


M7.5* 
MS'- 


2800 
2600 


0.014 
0.0035 


0.05* 
0.03 


na.b 

10" 


7±1 
7±1 


6.7±1.8 


<ir 

<8 


5.6+1 


11.9±1.7 

7.5±1 


11.9±0.7 


11.8±1.5 



" Jayawardhana et al. (2003) 

'' Comeron, Neuhauser & Kaas (2000) (A, = 0.11) 

■= Gizis(2002) 

'' upper limit refer to a 5cr detection threshold 

■^ smearing assumed (FWHM=larcsec) 



circum-sub-stellar material, as it appears only marginally in- 
dicative (for ChaHal) or even fails (for 2M1207). 

3. Disk Models 

In this section we explain the main features of the SED using 
the new TReCS photometry (combined with ISO photometry, 
if available) with models of BD disks. The models are based 
on the Chiang & Goldreich (1997) approximations, with mod- 
ifications from Dullemond, Dominik & Natta (2001) for flux 
conservation. The general methodology and application to BD 
disks is explained in Pascucci et al. (2003). In view of the un- 
certainties of the exact shape of the SED over a wider wave- 
length range due to the knowledge of only very few data points, 
we do not attempt to fine-tune the disk models with respect to 
variations of input parameters like disk inner and outer radii, or 
inclination. We merely want to concentrate the discussion on 
those parameters that have a most notable impact on the shape 
of the SED around lOjum. These are mainly the assumed disk 
geometry, and the dominant size and composition of the dust 
particle distribution. In the following, we assume standard as- 
tronomical silicates with opacities from Draine & Lee (1984). 
As we discuss only passive, re-processing disks, the assumed 
parameters of the BD mass and its luminosity are important. 
We directly use, or derive them from published values as indi- 
cated. 

3.1. ChaHal 

According to Comeron, Neuhauser & Kaas (2000), the spectral 
type of ChaHo-l is M7.5, and they estimate Te/f - 2770K and 
L = 0.01 ILo- Using the dusty BD model grids from Allard et 
al. (2000, 2001), we find that the available near-IR photome- 
try for this object is in much better agreement with an atmo- 
spheric model for a « IMyr old BD with T,ff = 2800K and 
L = 0.014Lq. We prefer these values for the BD photospheric 
parameters, which are anyway consistent with the Comeron et 
al. estimates, given their large error bars. Fig. |3 summarizes 
the results of our disk modeling in the MIR regime. The up- 
per panel refers to opacities of astronomical (amorphous) sili- 
cates of O.lfim size, while the lower panel displays the results 
of large, 2jjm sized grains of the same composition. Full lines 
indicate the SED expected for a flared model, in which the ver- 
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10 
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Fig. 3. SED and comparison of disk models for ChaHal. 
Measurements are diamonds, with associated error bars. The 
upper panels refers to small dust sizes (0.1/zm), the lower 
panel to larger grains (2/im). Flared (full line) and flat (dashed 
line) model predictions are compared. The contribution of the 
BD photosphere is approximated by a black-body model. The 
dashed-dotted lines refers to the disk model by Walker et al. 
(2004). Our T-ReCS measurements are shown with diamonds, 
while the ISOCAM data are indicated with boxes and the L' 
point with a triangle, respectively. 



tical scale height is derived from local hydrostatical equilib- 
rium at each disk radius. The dashed lines refer to the SED 
in a flat geometrical configuration. Contributions of the pho- 
tosphere are indicated as dotted lines. The dashed-dotted line 
indicates the results from Walker et al. (2004). We will refer to 
it later. 

We note that the L' photometry from Jayawardhana et al. 
(2003) aheady indicates a significant excess, and favors - ac- 
cording to our modeling - a flattish disk geometry. However, the 
superb quality of our new TReCS measurements, especially at 
I0.4fj.m, constrains the disk configuration to a standard, flared 
disk model. Standard, small, ISM grains can explain the sili- 
cate emission feature very well, although larger grains cannot 
be excluded. 
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Fig. 4. SED and comparison of disk models for 2M1207. 

3.2. 2MASS1 207-3932 

Following Gizis (2002), the spectral type of 2M1207 is M8, 
which translates to T^ff - 2600K according to the calibration 
of Luhman (1999). If we assume an age of lOMyr (see below), 
and scale the available JHK photometry to a distance of 70pc, 
we find an approximate bolometric luminosity of O.OOSSLq and 
a mass of O.OOSMq for this BD from the evolutionary tracks 
of Allard et al. (2000, 2001). These values are consistent with 
those estimated in MJBN. 

The interpretation of the SED for this BD is not unique. 
ISO observations are not available, and the two TReCS mea- 
surements, together with the upper limit at 7.7yum and the L' 
measurement clearly exclude a standard flared disk model, with 
ISM grain sizes (O.l/^m). Either a flat disk, and/or somewhat 
bigger grains (2/im) are compatible with the photometry. But 
large grains (up to 5;um) only seem to be possible in a flared 
configuration viewed from a high inclination angle (50deg). 
The flat geometrical configuration is, on the other hand, more 
robust with respect to inclination angle variations, because 
shadowing effects are less pronounced. 

Recently, Chauvin et al. (2004) report the detection of a 
giant planet candidate close to 2M1207. From near-IR spec- 
troscopy they derive a spectral type of L5-L9.5, and a likely 
temperature of T^ff - 1250±200Kfor this object, only 0.8 arc- 
sec separated from the primary. Although it should be possible 
to spatially resolve it with a 8m class telescope, we do not have 
any evidence for this object in any of our images (see Fig.|2}. 
But we note that the expected photospheric flux in the N-band 
is only around 0.2mJy, assuming a simple black-body scaled 
to their A'-band measurement. We therefore conclude that the 
excess flux we see is originating only from the disk material 
around 2M1207, and is not contaminated by the putative com- 
panion. 

4. Discussion 

Even with the new and improved MIR photometry an unam- 
biguous interpretation remains difficult, because grain sizes and 



dust spatial distribution are coupled. For both objects, flared 
models can in principle explain the observed MIR fluxes. A 
flared disk containing small grains heated in the upper surface 
layers is the most likely explanation for the pronounced sili- 
cate feature in ChaHal. A flat geometry, and substantial grain- 
growth are excluded for ChaHal. In the case of 2M1207 we 
find that a flared disk requires a large inclination and grains 
with sizes of a few //m to fit the two lOyum points. Equally wefl, 
a flat disk composed of small, or moderately large grains up to 
a: 2jum can account for the observations. 

We will now try to put these observations in an evolutionary 
scenario. 



4.1. Grain Growth and Age Evolution 

TWA membership of 2M1207 has been deducted both from 
proper motion (Gizis 2002), and radial velocities (MJBN). The 
age of the association itself is fairly well constrained: 8.3Myrs 
from dynamical age determination (Makaraov & Fabricius 
2001), ~ lOMyrs from isochrone fitting of its stellar popula- 
tion (Webb et al. 1999), and a similar age from Li abundance 
determination for late-type members (e.g. Soderblom et al. 
1998). However, the presence of optically thick circum-stellar 
disks in many of the prominent stellar members of TWA (TW- 
Hya, HD98800B, Hen3-600A, HR4796; see Jayawardhana et 
al. 1999) and signatures of active accretion like in the clas- 
sical TTS TW-Hya (Rucinski & Krautter 1983; Muzerolle et 
al. 2000) has always challenged a consistent evolutionary pic- 
ture for this group of stars. Gap formation and significant 
dust evolution are observed in all four stellar systems men- 
tioned above, and the detailed analysis of low-resolution N- 
band spectra allows to decompose the lOjum feature in distinct 
dust species: TW-Hya exhibits broad silicate emission, indica- 
tive for amorphous, large (2yum) olivine grains (Weinberger et 
al. 2002). HD98800B has highly processed dust, dominated 
by large amorphous olivine and crystalline forsterite (Schiitz, 
Meeus & Sterzik 2004). And also the spectrum of Hen3-600A 
features a rich mixture of crystalline silicate components, and 
even SiOa (Honda et al. 2003). 

In other words, all circum-stellar disks in TWA show signs 
of dust processing, and are in particular incompatible with 
a primordial dust size distribution like the ISM. Hence, the 
dust observed in these disks must have evolved on a lOMyr 
timescale since star formation. If the disk evolution is quali- 
tatively similar around stars and sub-stellar objects, then one 
might expect similar dust properties also in circum-sub-stellar 
disks. The indication of larger grains in 2M1207 from our mod- 
eling is therefore fully consistent with this notion. 

The case of ChaHal fits in this evolutionary scenario. Not 
much dust processing and/or growth has occurred, and one rea- 
son for this could well be its relatively young age (2 Myrs, 
or even < 0.1 Myrs, Gomez & Mardones 2003). Small grains 
are well-mixed throughout the vertical structure, which is it- 
self in hydrostatical equilibrium, and illuminated by the central 
source. 
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4.2. Comparison with ottier BDs 

Beside the two BDs discussed in this paper, the relation of disk 
geometry and grain size distribution has been investigated only 
in very few other BDs based on (ground-based) MIR photom- 
etry that sample the silicate feature well enough. 

Mohanty et al. (2004) discuss three BDs in the young 
(< IMyr) p-Ophiuchus star-forming region. For two of these 
objects (GY 310 and GY 11) they favor a standard, flared disk 
geometry, while a flat geometry cannot be ruled out for their 
third object (GY 5). Their data do not allow to draw strong 
constraints on the grain sizes, but small grains appear to bet- 
ter fit the silicate feature - at least in GY 310. In summary, 
these three, very young BD, do not show evidence for signif- 
icant grain and disk evolution, consistent with the notion that 
grain growth occurs at later stages, during disk evolution. 

But this simple, evolutionary, picture might be more diffi- 
cult to reconcile with the finding that other brown dwarf disks 
of similar age contain non-primordial dust, and larger grains 
then in the ISM. Good evidence for grain-growth is, e.g., pre- 
sented in Apai et al. (2004) for the ^IMyr old brown dwarf 
BD Tau 4. Although a precise age determination always suffers 
from notoriously uncertain assumption about sub-stellar evolu- 
tionary tracks and unknown distances, a much older age seems 
to be excluded because the object is already over-luminous for 
that particular isochrone (Martin et al. 2001). Nevertheless, this 
BD disk deviates from an equilibrated flared disk geometry, 
and rather suggests that a significant amount of grain growth 
and dust settling has already occurred in its short lifetime. Also 
in the case of ChaHQ;2, likely of similar age as ChaHcl, only a 
flat disk geometry (Apai et al. 2002), and/or a somewhat larger 
grain size distribution (Walker et al. 2004) can explain the ob- 
served SED. Both possibilities are indicative for some dust/disk 
evolution that has not - yet - happened in ChaHal . 

When, and how, does grain growth in circum(-sub)-stellar 
disks occur? 



4.3. Primordial Grain Growtii 

There is some observational evidence that larger grains (and 
therefore opacities having a shallow wavelength dependency) 
are observable very early in disks around T Tauri stars (see 
e.g. Wood et al. 2002). It has been suggested that grain growth 
might therefore happen already during the collapse phase, or 
very early in the formation phase of the circum-stellar disk 
(see, e.g. Suttner & Yorke 2001). Based on this assumption. 
Walker et al. (2004) model the SED of BD disks, too, and 
find reasonable agreement with those 12 BD disks that have 
ISOCAM MIR photometry. But not all of their model predic- 
tions are compatible with our observations. The Walker et al. 
(2004) model (indicated by a dot-dashed line in Fig. |3} fails, 
e.g., for ChaHcKl. Our new TReCS measurements rule out the 
possibility that the observed opacities (and therefore the grain 
sizes) deviate significantly from ISM properties. This means 
that at least for this objects grain growth, and dust processing, 
is not (yet) dominating the dust grain population as observed in 
flie MIR. 



4.4. Tine Role of Multiplicity 

The co-existence of dust disks with vastly different grain prop- 
erties has also been found in samples of classical TTS disks 
(Pryzgodda et al. 2003). Meeus et al. 2003, in particular, anal- 
yse the dust properties in a co-eval sample of TTS of simi- 
lar (late) spectral type in Chameleon. A detailed mineralogi- 
cal analysis of the N-band spectra of each individual object re- 
quires a different dust composition and size distribution, rang- 
ing from primordial ISM to highly processed dust. They note 
that the object that shows the highest degree of dust process- 
ing is in a close binary star, while the apparently unprocessed 
dust disk belongs to a single star. Without drawing strong con- 
clusions they speculate that the presence of a binary might in- 
fluence and help dust processing in disks. We might witness a 
similar trend in the observed BD disks in Chameleon. ChaHa2 
has been reported to be a close (0.2") binary (Neuhaeuser et al. 
2002), and its disk exhibits signs of significant dust processing 
and grain growth, in contrast to ChaHal, apparently a single 
BD (Joergens & Guenther 2001). 

There is no detailed model that explains the influence of a 
binary potential on the dust disk dynamics. But it seems plau- 
sible that a perturbing force acting on a circum-sub-stellar disk 
can induce enhanced collision, and therefore increase growth 
rates of dust particles in its disk (see, e.g. Dubrulle, Morfill & 
Sterzik 1995). 

In this light, the overluminosity of BD Tau 4 in the HR 
diagram, together with its dust properties, might be a conse- 
quence of a close companion. Interestingly, Itoh, Tamura & 
Nakajima (1999) find a faint companion candidate around this 
object, though at rather wide separation (^ 4.2"), and its phys- 
ical relation to BD Tau 4 is not established. 



5. Summary 

We report new high-quality MIR photometry from TReCS of 
two bona-fide brown dwarfs, and find excess emission in both 
cases, indicative for optical thick dust disks. We have inves- 
tigated the influence of diverse disk geometries (flared versus 
flat) and of different grain size populations (0. lyum versus 2yum) 
on the MIR-SED of these brown dwarfs. In particular, we find: 

- The MIR-SED of the ~ 2Myr old ChaHal is fully con- 
sistent with a standard flared disk model, and primordial 
(ISM) grain size properties. 

- The MIR-SED of the ^ lOMyr old 2M1207 can be better 
explained by the presence of larger dust, and/or a flat disk 
geometry. 

- The evolutionary timescale of dust in brown dwarfs appears 
to be similar as in TTS. 

- Similar to TTS, rapid grain growth and dust processing 
early in the BD disk evolution might be related to the pres- 
ence of a close companion. 

Acknowledgements. We thank C. Walker for providing her model 
fluxes of ChaHcfl in electronic form. Based on observations ob- 
tained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., un- 
der a cooperative agreement with the NSF on behalf of the 



Michael F. Sterzik et al.: Evolution of Young Brown Dwarf Disks in the Mid-Infrared 



Gemini partnership: the National Science Foundation (United States), 
the Particle Physics and Astronomy Research Council (United 
Kingdom), the National Research Council (Canada), CONICYT 
(Chile), the Australian Research Council (Australia), CNPq (Brazil), 
and CONICET (Argentina). We are indept to the TRecS/Gemini teams 
providing us excellent data from service observations." 



References 

Allard, F., Hauschildt, PH., Alexander, D.R., Tamanai, A., & 

Ferguson, J.W., 2000, in ASP Conf. Ser. 212, From Giant Planets 

to Cool Stars, ed. C.A. Griffith & M.S. Marley (San Francisco: 

ASP), 127 
Allard, F., Hauschildt, PH., Alexander, D.R., Tamanai, A., & 

Schweitzer, A., 2001, ApJ556, 357 
Apai, D., Pascucci, I., Henning, T, et al. 2002, ApJ, 573, LI 15 
Apai, D., Pascucci, I., Sterzik, M.F., et al., 2004, A&A, in press. 
Bontemps, S., Andre, P, Kaas, A.A., et al., 2001, A&A, 359, 269 
Bouwman, J., Meeus, G., de Koter, A., et al., 2001, A&A, 375, 950 
Chauvin, G., Lagrage, A.-M., Dumas., C, et al., 2004, A&A, in press 
Chiang, E.I., & Goldreich, P, 1997, ApJ, 490, 368 
Chiang, E.I., Joung, M.K., Creech-Eakman, M.J., et al., 2001, ApJ, 

547, 1077 
Comeron, R, Neuhauser, R., & Kaas, A.A., 2000, A&A, 359, 269 
Draine, B.T., & Lee, H.M., 1984, ApJ, 285, 89 
Dubrulle, B., Morfill, G., & Sterzik, M., 1995, Icarus, 114, 237 
Dullemond, C.P, Dominik, C, & Natta, A., 2001, ApJ, 560, 957 
Gomez, M., & Mardones, D., 2003, AJ, 125, 2134 
Gizis, I.E., 2002, ApJ, 575, 484 

Honda, M., Kataza, H., Okamoto, Y.K., et al., 2003, ApJ, 585, L59 
Itoh, Y, Tamura, M., & Nakajima, T., 1999, AJ, 117, 1471 
Jayawardhana, R., Hartmann, L., Fazio, G., et al., 1999, ApJ, 521, 

L129 
Jayawardhana, R., Ardila, D.R., Stelzer, B., & Haisch, K.E.Jn, 2003, 

AJ, 126, 1515 
Joergens, V., & Guenther, E., 2001, A&A, 379, L9 
Klein, R., Apai, D., Pascucci, I., Henning, T, & Waters, L.B.F.M., 

2003, ApJ, 593, L57 
Liu, M.C., Najita, J., & Tokunaga, A.T., 2003, ApJ, 585, 372 
Luhman, K.L., 1999, ApJ, 525, L466 
Makarov, V.V., & Fabricius, C, 2001, A&A, 368, 866 
Martin, E., Dougados, C, Magnier E., et al., 2001, ApJ, 561, L195 
Meeus, G., Waters, L.B.F.M., Bouwman, J., v.d.Ancker, M.E., 

Waelkens, C, & Malfait, K., 2001, A&A, 365, 476 
Meeus, G., Sterzik, M.R, Bouwman, J., & Natta, A., 2003, A&A, 409, 

L25 
Mohanty, S., Jayawardhana, R., Natta, A., et al., 2004, A&A, in press. 
Mohanty, S., Jayawardhana, R., & Barrado y Navascues, D., 2003, 

ApJ, 593, L109 (MJBN) 
Muench, A.A., Alves, J., Lada, C.J., & Lada, E.A., 2001, ApJ, 558, 

L51 
Muzerolle, J., Calvet, N., Briceno, C, Hartmann, L., & Hillenbrand, 

L., 2000, ApJ, 535, L47 
Neuhauser, R., Brandner, W., Alves, J., Joergens, V., & Cameron, R, 

2002, A&A, 384, 999 
Natta, A. & Testi, L., 2001, A&A, 376, L22 
Natta, A., Testi, L., Comeron, R, et al., 2002, A&A, 393, 597 
Pascucci, I., Apai, D., Henning, T., & Dullemond, C.P, 2003, ApJ, 

590, LI 11 
Persi, P, Marenzi, A.R., Olofsson, G., et al., 2000, A&A, 357, 219 
Pryzgodda, R, van Boekel, R., Abraham, P, et al,. 2003, A&A, 412, 

L43 
Rucinski, S.M., & Krautter, J., 1983, A&A, 121, 217 



Schutz, O., Meeus, G., & Sterzik, M.R, 2004, A&A, in press. 

Suttner, G., & Yorke, H.W., 2001, ApJ, 551, 461 

Telesco, C. M., Pina, R. K., Hanna, K. T, Juhan, J. A., Hon, D. B., 

Kisko, T M., 1998, SPIE 3354, 534 
van Boekel, R., Waters, L.B.F.M., Dominik, C, et al. 2003, A&A, 

400, L21 
Walker, C, Wood., K., Lada, C.J., Robitaille, T, Bjorkman, J., & 

Whitney, B., 2004, MNRAS, in press. 
Webb, R.A., Zuckerman, B., Platais, I., et al., 1999, ApJ, 512, L63 
Weinberger, A.J., Becklin, E.E., Schneider, G., et al., 2002, ApJ, 566, 

409 
Wood, K., Wolff, M.J., Bjorkman, I.E., & Whitney, B., 2002, ApJ, 

564, 887 



